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ABSTRACT: Future Higgs factories can reach impressive precision on Higgs property measure- 
ments. In this paper, instead of conventional focus of Higgs precision in certain interaction 
bases, we explored its sensitivity to new physics models at the electron-positron colliders. In 
particular, we studied two categories of new physics models, Standard Model (SM) with a 
real scalar singlet extension, and Two Higgs Double Model (2HDM) as examples of weakly- 
interacting models, Minimal Composite Higgs Model (MCHM) and three typical patterns 
of the more general operator counting for strong interacting models as examples of strong 
dynamics. We performed a global fit to various Higgs search channels to obtain the 95% C.L. 
constraints on the model parameter space. In the SM with a singlet extension, we obtained 
the limits on the singlet-doublet mixing angle sin 0, as well as the more general Wilson co- 
efficients of the induced higher dimensional operators. In the 2HDM, we analyzed tree level 
effects in tan 8 vs. cos(3—a) plane, as well as the one-loop contributions from the heavy Higgs 
bosons in the alignment limit to obtain the constraints on heavy Higgs masses for different 
types of 2HDM. In strong dynamics models, we obtained lower limits on the strong dynamics 
scale. In addition, once deviations of Higgs couplings are observed, they can be used to dis- 
tinguish different models. We also compared the sensitivity of various future Higgs factories, 
namely Circular Electron Positron Collider (CEPC), Future Circular Collider (FCC)-ee and 
International Linear Collider (ILC). 
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1 Introduction 


The discovery of a Standard Model (SM)-like Higgs boson at the Large Hadron Collider (LHC) 
strongly motivates the studies of its properties. Lepton colliders running at center of mass 
energies of around 240 GeV or above are ideal Higgs factory machines for the precision studies 
of the Higgs boson properties. Several plans for future lepton colliders, in particular, Higgs 
factories, have been proposed, including the Circular Electron Positron Collider (CEPC) in 
China [1], the electron-positron stage of the Future Circular Collider (FCC-ee) at CERN 
(previously known as TLEP [2]), and the International Linear Collider (ILC) in Japan [3]. 


te —> hZ, Higgs cross sections can be measured to 


Through the Higgsstrahlung process, e 
an astonishing precision: about 0.2%—0.5% for “inclusive” ete~ — hZ cross section and 
ete > hZ,h — bb exclusive cross section, and about a few percent for many other exclusive 


cross sections, under typical running scenarios. The Compact Linear Collider (CLIC), while 


only planning to run at energies 350 GeV and above, could also probe the Higgs couplings very 
well by measuring the WW fusion process at high energies [4, 5], as well as measuring the top 
Yukawa and the triple Higgs couplings. These colliders could provide exciting opportunities 
in Higgs studies and greatly improve our understanding of the physics at the TeV scale. 
Any deviation of the Higgs couplings from their SM predictions is a strong indication of new 
physics beyond the SM. Even if no observation is observed, it will provide tight constraints 
on new physics models. 

In the study of Higgs boson properties and the sensitivity of Higgs precision measure- 
ments to new physics models, two “model-independent” approaches are usually taken. One is 
the so-called “kappa” framework, in which «; is defined as the Higgs coupling to SM particles, 
normalized to its SM value. The kappa framework features the simple and direct connection 
to measurable cross sections and a straight forward inclusion of new light degree of freedom 
beyond the SM, e.g., exotic Higgs decays [6-8]. Deviation of k; from their SM expectation 
indicates contributions from new physics [9-20]. Another approach using the language of the 
Standard Model Effective Field Theory (SMEFT) for Higgs couplings has also been exten- 
sively studied, both at the LHC and at future colliders [21-36]. The fitting parameters are 
the coefficients of various higher dimensional operators. For both approaches, global fits to 
the SM observed values of ø x BR for various Higgs search channels are usually performed 
to obtain robust constraints and correlations on «k; couplings in the kappa framework, or 
Wilson coefficients in EFT framework. These fitting results are then further used to impose 
constraints on parameter spaces of a specific model when model parameters are mapped to 
either those k; couplings or Wilson coefficients. 

Such model independent approaches are very useful in evaluating the performance of 
measurements and providing general constraints on new physics models. However, they have 
certain disadvantages. The precision reaches of the global fits in such approaches usually 
suffer from the large level of degeneracy among fitting parameters. For specific models, 
the parameter space is usually much smaller, with various Higgs couplings determined by 
a small set of parameters. Furthermore, while the coupling fit results are provided in the 
official documents for each of the colliders, the corresponding correlation matrix is often not 
provided. For poorly determined quantities, linear approximation might not be sufficient. 
Directly applying these coupling fitting results alone on specific models without including 
correlation matrix may result in overly conservative estimations. A specific example of this 
is shown later in Section 5. Moreover, in most model independent approaches, certain model 
assumptions and simplifications are usually made, which may not be valid for a particular new 
physics model. Typical examples are the omission of Higgs couplings with Lorentz structures 
different from the SM ones (e.g. hZ,,Z"” vs. the SM hZ"Z,,) in the kappa framework, or 
the absence of possible light beyond the SM (BSM) states in SMEFT. The latter may also 
lead to Higgs exotic decays, which is hard to study in a model independent way, while many 
appealing BSM scenarios with light states can be sufficiently probed at e~e* colliders [7]. 
These disadvantages severely limit the usefulness of model independent methods when applied 
to specific BSM models. 


To demonstrate the sensitivity of Higgs factories to BSM new physics models and access 
the capability of various machine options, it is thus important to directly study the implication 
of Higgs precision measurements on model parameter space. We focused on a few generic 
types of new physics models in our study. For weakly coupled models, we used the SM plus 
a real scalar singlet, and Two Higgs doublet Model (2HDM) as two prototypes, which are 
incorporated as the Higgs sector for many BSM scenarios such as supersymmetric models, 
and left-right symmetric models. In particular, we studied the tree level effects in general 
case, as well as the one-loop contributions from the heavy Higgs bosons. For models with 
strong dynamics, we studied the Minimal Composite Higgs Model (MCHM), including all ten 
different Fermion embeddings, as well as three typical patterns of the more generic operator 
counting for strong interacting models with a light Higgs. With these two classes of Higgs 
extension models covering both weakly interacting theories and strong dynamics, we could 
develop a better understanding of the physics potential of future Higgs factories. 


We performed a global fit to various Higgs search channels in those two classes of specific 
models. In the case of no deviations observed in the future Higgs factories, a 95% C.L. 
constrained parameter space was obtained. For the SM with a singlet extension, we obtained 
limits on the singlet-doublet mixing angle sin 0, as well as the more general Wilson coefficients 
of the induced higher dimensional operators. For the 2HDM tree level effects, it tells how much 
deviation from the alignment limit we can tolerate. For loop effects under the alignment limit, 
a constrained range of heavy Higgs masses can be obtained. In the MCHM, we translated 
the Higgs precision measurements into constraints on the vacuum misalignment parameter 
€, which further constrain the composite scale f. In the three cases of strong interacting 
models with a light Higgs, we obtained the lower limits on the strong dynamics scale m, as 
a function of the strong interaction coupling gą. In cases of an observed deviation of the SM 
Higgs couplings, our studies also demonstrated how to distinguish different models. 


While our analyses mostly used the CEPC results, comparison with the reach of FCC- 
ee and ILC was also performed. Run-I, Run-II and high luminosity LHC Higgs precision 
measurements were included as well. We also compared the sensitivity of fitting using A(o x 
BR)/(o x BR)gm of variously Higgs search channels directly, with the fitting results using the 
constraints on effective couplings «; from existing studies. 


The rest of this paper is organized as follows. In Section 2 we summarized the run sce- 
narios and the estimated precisions of Higgs measurements for various future lepton collider 
Higgs factories, gathered from the corresponding official documents. These inputs were used 
to obtain constraints in the model parameter spaces. In Section 3, we presented the global 
fitting method used in our analyses. In Section 4, we started with the simple case of the 
SM plus a singlet as a warm up. The results for Two Higgs doublet Models were presented 
in Section 5 and the results for Composite Higgs Models were presented in Section 6. We 
concluded in Section 7. The LHC Run-I Higgs measurement results and the projected preci- 
sions for future LHC runs were collected in Appendix A. We also listed the formulae of loop 
corrections to various Higgs couplings in the 2HDM along the alignment limit in Appendix B. 


2 The Higgs precision measurements at future lepton colliders 


At future lepton colliders, the dominant channel to measure the Higgs boson properties is the 
Higesstrahlung process, ete — hZ, at center of mass energies of around 240—250 GeV. Due 
to the nature of lepton colliders, both the inclusive cross section, o (hZ), and the exclusive ones 
of different Higgs decays in terms of o(hZ) x BR, can be measured to remarkable precisions. 
The invisible decay width of the Higgs can also be very well constrained. In addition, the 
cross section of WW fusion process for Higgs production grows with energy. While it can not 
be measured very well at 240—250 GeV, at higher center of mass energies (in particular at 
linear colliders), such fusion process becomes significantly more important and can provide 
crucial complementary information. For ys > 500 GeV, tth production can also be used as 


well. 
collider CEPC | FCC-ee ILC 
V/s 240 GeV | 240GeV | 250 GeV 350 GeV 500 GeV 
f Lat 5ab-' | 10ab-1 | 2 ab"? 200 fb~! 4 abo? 
production Zh Zh Zh Zh voh Zh voh tth 
Ac/o 0.51% 0.4% 0.71% | 2.1% - 1.06 - - 
decay A(o- BR)/(o- BR) 
h —» bb 0.28% 0.2% 0.42% | 1.67% | 1.67% | 0.64% | 0.25% | 9.9% 
h > c 2.2% 1.2% 2.9% | 12.7% | 16.7% | 4.5% | 2.2% - 
h —> gg 1.6% 1.4% 2.5% | 9.4% | 11.0% | 3.9%] 1.5% : 
h —> WW~* 1.5% 0.9% 1.1% | 8.7% 6.4% | 3.3% | 0.85% - 
h> Tr 1.2% 0.7% 2.3% | 4.5% | 244% | 1.9% | 3.2% - 
h-> Z2Z* 4.3% 3.1% 6.7% | 28.3% | 21.8% | 8.8% | 2.9% - 
h> yy 9.0% 3.0% | 12.0% | 43.7% | 50.1% | 12.0% | 6.7% = 
h > pty 17% 13% 25.5% | 97.6% | 179.8% | 31.1% | 25.5% - 
(vi)h — bb 2.8% 2.2% 3.7% - - - - - 


Table 1. Estimated statistical precisions for Higgs measurements obtained at the proposed CEPC 
program with 5 ab~! integrated luminosity [1], FCC-ee program with 10 ab~! integrated luminosity [2], 
and ILC with various center of mass energies [37]. 


To set up the baseline of our study, we hereby listed the run scenarios of various machines 
in terms of center of mass energy and the corresponding integrated luminosity, as well as the 
estimated precisions of relevant Higgs measurements that we used in our global analyses: 


e CEPC According to the preCDR [1], CEPC plans to collect 5ab~! data at 240 GeV. 
The estimations on the measurements of the Higgsstrahlung process ete~ — hZ with 
various final states, as well as the WW fusion process with Higgs decaying to bottom 
pairs (ete~ —> vih, h — bb), are summarized in Table 1. 


e FCC-ee The FCC-ee CDR is expected to be finished by 2018 [38]. At the current 
moment, the TLEP whitepaper [2] still provides the most updated estimations, which 
assumes total luminosities of 10ab~! at 240 GeV and 2.6ab~! at 350 GeV. The esti- 
mations for the ete~ — hZ measurements at 240 GeV, as well as h — bb channel in 
WW fusion are listed in Table 1. Assuming statistical uncertainties only, Ref. [2] also 
estimated that the WW fusion process with h — bb can be measured to a precision of 
0.6% at 350 GeV. We did not include 350 GeV WW fusion process in our global fit 
since it has little impact on the constrained model parameter spaces considered in our 
study. 


e ILC The proposed run scenarios in the ILC TDR [3] have been updated in re- 
cent documents [37, 39], which suggested that the ILC could collect 2ab~! data at 
250 GeV, 200fb~! at 350GeV, and 4ab7! at 500GeV. However, the estimations of 
signal strengths, summarized in Ref. [37], are only available for smaller benchmark lu- 
minosities for which the full detector studies were performed. We took these estimations 
and scaled them up to the current run scenarios, assuming statistical uncertainties dom- 
inate [40]. The scaled estimations are summarized in Table 1. Such scaling provides a 
reasonable approximation as long as the luminosities are not excessively large and the 
systematic uncertainties are under control. 


With large center of mass energies up to 3 TeV, CLIC is also able to measure the Higgs prop- 
erties very well through the WW fusion process [4, 5]. On the other hand, with its extensive 
coverage of energy scales, the primary goal of CLIC is to directly search for new particles, in 
particular the ones coupled to SM particles only through electroweak interactions. CLIC is 
advantageous in the possibility of directly producing new particles that modifies Higgs prop- 
erties at low energy, and a large span of observables at different scale. Consequently, a study 
on the CLIC physics potential for various models would require additional considerations be- 
yond Higgs precision physics. A comprehensive study of the CLIC physics potential including 
both the direct and indirect searches of new physics is beyond the scope of this paper. 

In our global fit to the Higgs measurements, we only included the rate information for 
the Higgsstrahlung as well as the WW fusion process. The measurements of the angular 
distributions in the Higgsstrahlung process can provide important information in addition to 
the rate measurements alone [29, 30]. The diboson process, ete” + WW, can be measured 
to a great precision, which imposes very strong constraints on the anomalous triple gauge 
couplings (TGCs) [41-44]. These measurements are very helpful in probing new physics. 
In particular, it is shown in Ref. [34] that the inclusion of these measurements are crucial 
for constraining new physics in a global EFT framework, which exhibits large flat directions 
with the Higgs rate measurements alone. However, their impacts are significantly smaller 
for specific models like the SM plus a singlet, 2HDM and MCHM, due to a much smaller 
model parameter space (compared with SMEFT). Therefore, they were included in our global 
fit of operator approach of strong dynamics models only. The electroweak (EW) precision 
measurements at the Z-pole also impose strong constraints on the new physics [45, 46]. The 


current constraints from the Large Electron Positron Collider (LEP) could be significantly 
improved by a Z-pole run at any of the future lepton colliders. While these constraints were 
not explicitly considered in our study, we did restrict ourselves to models with suppressed 
EW precision corrections (e.g. by imposing custodial symmetries) such that these constraints 
are automatically satisfied. 

It is also important to compare the reaches of the future Higgs factories to that of the 
LHC. For the LHC Run-I Higgs measurements with 5 fb~! integrated luminosity at \/s = 7 
TeV and 20 fb—! at \/s = 8 TeV, we used the results in Ref. [47]. For the LHC with 300 fb! 
and 3000 fb~! luminosities, we used the ATLAS projections in Ref. [48], which collects the 
information from several other studies. The detailed inputs are listed in Appendix A, with 
the LHC Run-I results in Table 8 and the ATLAS projections for LHC 300 fb~! and 3000 fb~! 
summarized in Table 9. 


3 Global fit framework 


To transfer the estimated error on the experimental measurements to the constraints on the 
model parameters, we made a global fit by constructing the x? with the profile likelihood 
method 


py ey en 
o2, 
a Hi 
Here pBSM — oan for various Higgs search channels and o,,, is the experimental pre- 


cision on a particular channel. We note that the correlations among different o x BR are 


usually not provided, and are thus assumed to be zero in the fits. PSM 


i is predicted in each 


specific model, depending on model parameters. For the LHC Run-I, the measured er? and 
obs 


corresponding op; are given in Table 8. In our analyses, for the future colliders, u; 


; are set 


to be the SM value: > = 1, assuming no deviation to the SM observables are observed. 
The corresponding o,,, are the estimated error for each process, as shown in Table 1 for the 
CEPC, FCC-ee, ILC and Table 9 for the LHC. For the ILC with three different center of 
mass energies, we summed the contribution from each individual channel. For one or two 
parameter fit, the corresponding Ay? = x? — x2,,, for 95% C.L. is 3.84 or 5.99, respectively. 

We fitted directly to the signal strength ui, instead of the effective couplings «xi. The 
latter are usually presented in most experimental papers. While using «-framework is easy to 
map to specific models, unlike u;, various k; are not independent experimental observables. 
Ultimately, fitting to either u; or ki should give the same results, if the correlations between 
ki are properly included. Those correlation matrices, however, are typically not provided. 
Therefore, fit to «; only, assuming that they are uncorrelated, usually leads to a more relaxed 
constraints. Comparison of u-fit versus «-fit results is given later in the example of the 2HDM. 

In the last example of our study for future collider constraints on generic strong dynamics, 
we adopted the fitting results in EFT coefficients given in Ref. [34] since various scenarios 
can only be meaningfully discussed in power counting of the structure of the induced EFT 
operators. EFT provides more complex and rich structure for the Higgs couplings, which 


requires the inclusion of additional measurements to ensure the parameters are reasonably 
well-constrained. In particular, diboson process (ee — WW) and the angular observables 
in ete — hZ were included in the global fit in addition to the Higgs rate measurements. 
In the EFT framework, the y? from the Higgs rate measurements and other measurements 
can be constructed in a similar way as Eq. (3.1), with uP™ being a function of the EFT 
parameters (7.e., Wilson coefficients). 


4 SM with a real singlet extension 


We first applied this global fit to the simplest extension to the SM with a real scalar singlet. 
The general Lagrangian for this model is, 


Phage Las) z sms? — Ag S(HtH) — 5 AsHS?(HTH) : Ass’ 7 asst, (4.1) 
where H is the SM Higgs doublet and § is the new real singlet field. 

This simplest extension to the SM could already induce many test scenarios. Such model 
is of particular interest as it can in general help enhance the electroweak phase transition to 
be a strong first order one [49], which is needed for electroweak baryogengesis. With certain 
constraints on the model parameters, it is also a good description of the scalar section for the 
next to minimal supersymmetric Standard Model [50, 51] in the decoupling regime. 

The Lagrangian given in Eq. (4.1) can be categorized into two scenarios: Zə preserving 
and Zə breaking. In the Zə preserving scenario, the Lagrangian is further simplified: Asg = 
Ag = 0, reducing the number of free parameters in this model to be three. Note that in 
the Zə limit, it is still possible to have spontaneous Z2 breaking once S acquires a vacuum 
expectation value (VEV), leaving us same number of terms for the interaction Lagrangian 
comparing to the Z breaking case, but still with only three free parameters. Since the 
purpose of this study is to focus on the Higgs physics implications, instead of the extraction 
of singlet model parameters, we did not single out spontaneous Z breaking as a separate 
scenario. 

The biggest difference between the Zə preserving and the Z2 breaking is that the latter 
enables the singlet-Higgs mixing. The 125 GeV Higgs is one of the mass eigenstates, which 
is a mixture of real component of the SM doublet Asm and the singlet S: 


hia5 = cos hsm + sind S, (4.2) 
where 0 is the mixing angle. All of the SM-like Higgs couplings to other SM particles are 
scaled down by cos at tree level: ki = g? Mkeimglet / gM = cos 0. Since cos 0 does not exceed 


unity, the modifications of the Higgs couplings through mixing with the singlet always lead to 
a reduction of the Higgs couplings to SM particles. This mixing angle description remains the 
same for both a light and a heavy singlet-like scalar and is used widely due to its simplicity. 

As the first and simplest application of the global fit, we applied the global A? fit to the 
SM plus a singlet model with the only fitting parameter being sin 0. The Ay? distribution is 
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Figure 1. For SM plus one real singlet model, left panel shows the Ax? as a function of the singlet- 
SM Higgs doublet mixing angle |sin6| in the Z)-breaking case for CEPC (red), FCC-ee(green) and 
ILC (blue) from Higgs precision measurements. Right panel shows the 95% C.L. limit on the Wilson 
coefficient cgv?/A? vs. cyv?/A?. The vertical dashed lines indicate the one parameter fit limit on 
cyv’/A? with ce set to be 0. The corresponding values of the mixing angle |sin6| and changes to 
trilinear Higgs coupling A)3/ ASM are shown in the upper and right axes of this figure. 


shown in the left panel of Fig. 1 for the CEPC, ILC and FCC-ee precisions. The mixing angle 
| sin 0| is constrained to be 0.62, 0.058 and 0.052 for the CEPC, ILC and FCC-ee, respectively, 
at 95% C.L. 

While the mixing angle sin 0 captures the most important tree-level effect for the Higgs 
properties, it does not characterize the changes to the Higgs trilinear coupling, and neither 
the loop corrections to Higgs physics from the singlet field, e.g. Ref. [52]. To fully explore 
the changes of the SM Higgs property in models with an extra real scalar singlet, we adopted 
the EFT language to examine all the possible effects. 

After integrating out the singlet field, the general EFT with dimension-six operators can 
be written as, 

Ac= Hoy + S06, with On = 5 (Onl Ht H)?, Os = |HtH. (4.3) 
The operator Oy induces a universal shift in the Higgs couplings through the Higgs wave- 
function renormalization. The mixing angle 0 can be mapped with cy as 1 — cos 0 ~ 67/2 ~ 
1/2 x cyv?/A?, for v = 246 GeV. 

The Wilson coefficients for Oy and Og can be mapped through the tree-level processes 
in the Z breaking scenario as [53, 54] 


CH _ Ay C6 _ ( ASHA y Sst) (4.4) 


4 6 
2M% 66mg 


and through loop-level processes in both the Zə preserving and the Z breaking scenario 
as [54], 


cH ASH co ASH (4.5) 


The sign of the Oy operator is positive semi-definite, for both Z2 breaking and Z2 preserving 
scenario, which is consistent with the mixing angle description that cos@ < 1 and the Higgs 
couplings are always reduced comparing to that of the SM. The Og operator modifies mostly 
the Higgs couplings to electroweak gauge bosons and top quarks at loop-level. In our analysis 
we took into account the loop corrections to Higgs to weak gauge boson couplings via Og 
operators using calculations detailed in Ref. [55]. 

We performed our global fit to this singlet scenario with CEPC, ILC, and FCC-ee Higgs 
precision measurements, with 95% C.L. limits in the cgv?/A?-cyv2/A? plane shown in the 
right panel of Fig. 1.! The corresponding values of the mixing angle sin@ are shown in the 
upper axis of the figure, and the corresponding relative change to the SM trilinear Higgs 
coupling AA3/ AŞM is shown in the right axis of the figure. The one-parameter fit results with 
cyv/A? (or equivalently, sin @) are also shown in the vertical dashed lines. 

From Fig. 1, the allowed range in cgv?/A? is almost one order of magnitude larger than 
that of the cyv?/A?. In the general case, at 95% C.L., the maximum allowed Wilson coefficient 
cyv"/A? is 0.032, 0.021 and 0.005, corresponding to the max value of the singlet-double 
mixing parameter sin @ of 0.18, 0.15 and 0.07, for the CEPC, FCC-ee and ILC, respectively. 
The loop-modifications to Higgs to gauge boson couplings through the modifications of the 
trilinear Higgs coupling have sizable energy dependence. Hence, the ILC reaches a much 
better precision in this two-parameter fit via measuring the Higgs processes through different 
production modes well at separated center of mass energies of 250 GeV and 500 GeV. On the 
other hand, in the restrict one parameter fit case with cg set to zero, the limit on cgv?/A? 
improves to 0.0038, 0.0034 and 0.0028, for the CEPC, ILC and FCC-ee, respectively. These 
limits on cyv?/A? correspond to the maximum value of the singlet-double mixing parameter 
sin 0 derived from the left panel of this figure.” 


5 Two Higgs doublet model 


2HDMs are very generic BSM Higgs extension of the SM, including MSSM, gauge extensions 
(such as Left-Right symmetric model), and flavor models [58, 59]. Understanding the Higgs 
physics potential of 2HDM at future lepton colliders provides unique information covering a 
broad class of BSM. 


‘It should be noted that, the ILC 500 GeV could directly probe the Og operator by measuring the double 
Higgsstrahlung process, e*e~ — Zhh, which is not included in our analysis. See Refs. [36, 56] for more details. 

?The electroweak precision observables also provide relevant constraints on the parameter space of this 
model, mainly from the running and mixing of the Oy operator. For detailed analysis in various cases, see 


e.g., Refs. [53, 57]. 


Two SU(2)z scalar doublets ®;, i = 1,2 are introduced in 2HDM, 


oF 
i a +494 one) : (51) 


Each obtains a VEV v or ve after electroweak symmetry breaking (EWSB) with v? + v3 = 
v? = (246 GeV)?, and v1/v2 = tan £. 
The 2HDM lagrangian for Higgs sector can be written as 


£L=)°|D, i)? — V (81, 2) + Lyuk, (5.2) 
i 
with the Higgs potential 
V (1, 82) = m7, 610, + m3, 0462 — m (DiD + h.c.) + Maja) ae 2 aiaa)? 
HACO] D1) (D122) + M(E} D) (001) + 5 Doaa)? + he], 6.3) 
assuming CP-conserving and a soft Z2 symmetry breaking term mîz. 
After EWSB, one of the four neutral components and two of the four charged compo- 


nents are eaten by the SM Z, W~, providing their masses. The remaining physical mass 
eigenstates are the two CP-even Higgses h and H, with mp < my, one CP-odd Higgs A, 


as well as a pair of charged ones H=. Instead of the eight parameters appearing in the 
Higgs potential: (m? , M22, M?s, À1,2,3,4,5), a More convenient choice of the parameters is: 
(v, tan 8, a, Mh, MH, MA, Mg++, M?3), in which a is the rotation angle diagonalizing the CP- 


even Higgs mass matrix. 

The CP-even Higgs couplings to the SM gauge bosons are: gavy « sin(8 — a), and 
guvv xX cos(ĝ — a). For cos(8 — a) = 0, the light CP-even Higgs h couples to the gauge 
boson with full SM coupling strength while H decouples. This is usually referred to the 
“alignment limit” [60] with h identified as the SM Higgs °. For sin( — a) = 0, the opposite 
occurs with the heavy H being identified as the SM Higgs. While it is still a viable option 
for the heavy Higgs being the observed 125 GeV SM-like Higgs [61], the parameter spaces 
are squeezed with the tight direct and indirect experimental constraints. Therefore, in our 
analyses below, we identified the light CP-even Higgs h as the SM-like Higgs with mp, fixed 
to be 125 GeV. 

The most general Yukawa interactions of ®1 2 with the SM fermions under the Z2 sym- 
metry is 

— LYuk = YuQrio ur + YuQ,Padp + YD, ®.ep + h.c. (5.4) 


where ude are either ®; or ®2. Depending on the interactions of ®; coupling to the fermion 
sector, there are typically four types of 2HDM: 


e Type-I: ©; couples to all the fermions while ®2 does not couple to fermions at all. 


3Yukawa couplings of h to the SM fermions are also identical to the SM predictions under alignment limit. 
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e Type-II: ®; couples to the up-type quark, and ® couples to the down-type quarks and 
leptons. 


e Type-L: ®; couples to the quarks and ®2 couples to the leptons. 


e Type-F: ı couples to the up-type quarks and leptons while ®2 couples to the down-type 
quarks. 


For a review on different types of 2HDM as well as the phenomena, see Ref. [62]. 


Normalized Higgs couplings 
Kh Ki Ky, KH Ki KH KA A KA 
Typel | g9 gg mg dg Hg cotf -cotf -cot 
Typell | 9 -iig -g dig gg ge cot tanp tanp 
Type L| g9 æg dng sna dug eg cotf cot tan 
TypeF | gg -äng mg dag cg dag cotf tanp  — cotf 


Table 2. Higgs couplings to the SM fermions in the four different types of 2HDM, normalized to the 
corresponding SM value. 


After EWSB, the effective lagrangian for the light CP-even Higgs couplings to the SM 
particles can be parameterized as 


2m 
L= nn Zz, Zh + kw am we We he eye 


-Gi Gh + kyz Z Aw Al 


Br 
Bo Aw 2h — (ne > fft D CPAM y LA) , (5.5) 
f=u,c,t f=d,s,b fH=e,p,T 


in which «; is the SM-like Higgs coupling normalized to the corresponding SM value: ki = 
of2™ /g7¥ (same as x$ in Table 2). Table 2 summarizes all the tree-level non-zero «K; of 
the light CP-even Higgs h for the four different types of 2HDM, as well as the normalized 
couplings of the non-SM Higgses H and A. Note that while xg and k} are zero at tree level 
for both the SM and 2HDM, they are generated at the loop level. In the SM, both kg and x. 
receive contribution from fermions (mostly top quark) running in the loop, while xy receives 
contribution from W-loop in addition[54]. In the 2HDM, the corresponding hf f and hWW 
couplings that enter the loop corrections need to be modified to the corresponding 2HDM 
values. Expressions for the dependence of Kg, ky and Kz, on Ky and «s can be found in 
Ref. [63]. There are, in addition, loop corrections to Kg, Ky from extra Higgses. 

In the alignment limit of cos(8 — a) = 0, the light CP-even Higgs couplings are exactly 
identical to the SM ones. However, deviations of the Higgs couplings from the alignment limit 
are still allowed given the current LHC Higgs measurements [14, 64, 65]. Note that all the 
tree-level deviations to the SM couplings in Table 2 are parametrized by only two parameters: 
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@ and a, or more conveniently, tan 8 and cos(3 — qa). In sec. 5.1, we studied the constraints 
on the 2HDM tree level effects, namely, in the parameters space of tan 8 vs. cos(3 — a) with 
the Higgs precision measurements. 

There are, of course, loop corrections to all the «js above with non-SM heavy Higgses 
running in the loop. While their contributions are typically small comparing to the tree devi- 
ations, in the alignment limit, their contributions could be manifest given the high precision 
Higgs coupling measurements achieveable at future Higgs factories. The masses of the heavy 
Higgses will enter the loop corrections, as well as the soft Z2 symmetry breaking parameter 
m?,. We studied the constraints on the 2HDM loop effects in Sec. 5.2 under the alignment 
limit. 


5.1 2HDM tree-level results 


Performing a global fit to the Higgs rate measurements at the LHC as well as the CEPC*, we 
obtained the 95% C.L. region in the tan 8 vs. cos(3 — a) plane for various types of 2HDM, 
as shown in Fig. 2. For Type-I 2HDM with tan 6 = 2, | cos(3 — a)| is constrained to be less 
than about 0.5 with the current LHC Run-I data. With full LHC luminosity of 300 fb~', the 
range in |cos(G — a)| can be shrunk by about a factor of 2. At the HL-LHC with 10 times 
luminosity, the range can be further constrained to be less than 0.2. At the CEPC with 5 
ab™t, | cos(8—a)| is constrained to be less than about 0.08. The preferred range in cos(8 — a) 
is quickly shrunk at small tan 8, given that 
cos Q 


1 
AKu,d,e = sin l= 5 cos?(8 — a) 4 


The deviation of the Yukawa couplings from the experimental centre value is proportional to 


cos( — a) 


tanp (5.6) 


1/tan 8, resulting a reduced survival parameter space at small tan 8. The up-type Yukawa 
couplings are the same for all the four types of 2HDM. Therefore, the small tan 6 behavior is 
similar. 
For the Type-II, Type-L and Type-F 2HDM, the correction to the down quark and/or 
lepton Yukawa couplings receive tan 8 enhancement at large value, 
sin 8 


Ate = ET i= ; cos”(3 — a) — cos(8 — a) x tan B. (5.7) 


Here x is d,e in the Type-II, e in the Type-L and d in the Type-F. Therefore, the survival 
parameter space at large tan 8 is reduced significantly in all these three types. 

For the Type-II at the upper right panel of Fig. 2, as a result of larger tan 6 enhancement 
from Akg and small tan 8 enhancement from Aku, the region around tan 6 = 1 accommo- 
dates the largest deviation from the alignment. The current LHC Run-I constrain | cos(G—a)| 
less than 0.1 around tan 8 = 1 (except for the wrong-sign Yukawa couplings region [66, 67]), 
and the LHC Run-II 300 (3000) fb~* can reduce it to be less than 0.03(0.015). CEPC preci- 
sion measurements further reduce it to be less than 0.006. Similar behavior appears in Type-L 


“For future measurements, we have assumed that no deviation to the SM values is observed. 
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Figure 2. The allowed region in the plane of tan 8 vs. cos(8 — a) at 95% C.L. for the four types 
of 2HDM, given LHC and CEPC Higgs precision measurements. For future measurements, we have 
assumed that the measurements agree with SM predictions. The special “arm” regions for the Type-II, 
L and F are the wrong-sign Yukawa regions. See text for more details. 


and Type-F, with the small difference mostly coming from kẹ and «, parameter dependence. 
A summary of the 95% C.L. allowed maximum | cos(( — a)| range is given in Table 3. 


Because of those large or small tan 8 enhanced Higgs couplings deviation, we can examine 


—13- 


LHC Run-IT | LHC Run-Il 
Type LHC Run-I 300 fb-! 3000 fb-! CEPC 
Type-I tan 8 = 5 0.5 0.27 0.2 0.08 
Type-I] tan 6 ~ 1 0.1 0.03 0.015 0.006 
Type-L tan ~ 1 0.2 0.1 0.06 0.007 
Type-F tang ~ 1 0.1 0.03 0.02 0.005 


Table 3. Maximally allowed |cos(8 — a)| range at 95% C.L. given LHC and CEPC Higgs precision 
measurements. 


which Higgs decay channel provides the best constraints in certain region. Conservative 
7-parameters fit with CEPC Higgs measurements shows that oky < 0.16%, Ok, < 1.2%, 
One < 1.6%, on, < 1.3% [1]. 


Ky and k, typically provide the strongest constraints. For 
the Type-I, hbb coupling provides the strongest constrains at tan 8 < 1, while hZZ coupling 
constraints dominate for tan 8 = 5. For Type-II, hbb coupling dominates at large tan 8 and 
hgg coupling dominates at small tan 8. For Type-L, while the small tan 6 case is similar to 
that of Type-I, hrr dominates at large tan 8. Type-F is very similar to Type II, only that 
hrr enters at small tan 8 instead. 

The special “arm” region survived for the Type-I], L and Type-F in Fig. 2 is the so-called 
the wrong-sign Yukawa couplings [66], coming from fermion couplings with ®,: cosa/sin 2, 
as shown in Eq. (5.6). For tan 8 = 2/cos((G — a), there is no deviation of the corresponding 
couplings, which corresponds to exactly the “arm” central lines. For regions near this line, 
Ax flips its sign while the absolute value is small. The «k; for the other fermions that couple 
to ®2 is near —1, therefore “wrong-sign” while still surviving the Higgs measurements at 
the LHC. Actually, there would be another “arm” to make zero deviation for Eq. (5.7): 
tan 8 = -4 cos(3 — a). For Type-I, there is no such kind of coupling, and for other types, 
such region does not show up in our plot since we consider region tan 6 > 0.1. Smaller region 
of tan 8 is tightly constrained by the perturbativity of top Yukawa couplings. 

To compare the potential reach for different Higgs factories, we show the 95% C.L. reach 
in tan 6 vs. cos(G — a) plane based on the estimated Higgs measurement precision at CEPC 
(red), ILC (blue), and FCC-ee (green) in Fig. 3, for four different types of 2HDM. The red 
CEPC region is the zoom-in of the red region in Fig. 2. The 95% C.L. green and blue regions 
are almost the same, showing that ILC and FCC-ee has about the same constraining power, 
and both are slightly better than the CEPC results, with about 70% — 90% of the maximum 
| cos(8 — a)| range. 

In Fig. 3, we also studied the CEPC results with different luminosity to get a knowledge 
of various running scenarios. With the dashed red lines from the outer to inner, we presented 
the CEPC 2.5, 10, and 25 ab~! luminosity reach. In particular, results with CEPC 10 ab~! 
are almost the same as the FCC-ee 10 ab™t. 


We also show the comparison between the results using the signal strength p-fit (red 
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Figure 3. The comparison between the CEPC (red region), ILC (blue region) and FCC-ee (green 
region) reach in the plane cos(8 — a) vs. tan 8. A tiny arm region for Type-L is omitted for clarity. 
We also show the global fitting results to effective couplings from the 7 parameter fit of CEPC, instead 


of fitting to Au/p, in red solid line. Scaled CEPC results with 2.5 ab~!, 10 ab~', 25 ab! are shown 
in dashed lines, from outer to inner region. 


region) or the effective coupling «-fit (red line) for the CEPC precision. We adopted the 
precision for k using CEPC 7-parameter fit [1]. The results with «-fit are less restrained than 


that of the p-fit since no correlations between «; have been taken into account. Numerically, 
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except for the Type-I, the range of | cos(G — a)| at tan 8 ~ 1 with «-fit is about 3 times that 
of the y-fit for all three types of 2HDM. For Type II and Type-F, the CEPC «-fit results 
are similar to those of the HL-LHC, while for Type-I and -L, the CEPC «-fit results are 
still better than those of the HL-LHC. This demonstrates the under estimation of the Higgs 
physics potential if using « results without full correlation information. 


Model Distinction, cos(8-a)<0 Model Distinction, cos(6-a)>0 


0.00 


A(ky/kg)=(ky/kg—1) 


-0.05 


3, 
~0.19 50 0.05 0.00 0.05 0.10 0.10 -0.05 0.00 0.05 0.10 


Figure 4. A(ky/ka) vs Are for four different types of 2HDM, varying tan 8 and cos(8 — a). The left 
panel is for cos(8 — a) < 0 and the right panel is for cos(8 — a) > 0. The green, blue and red lines 
are for |cos(8 — a)| = 0.01, 0.02, 0.03 respectively. The dashed lines are for different tan 8 values, as 
labeled in the plot. The black cross indicates the estimated experimental errors with a random central 
point. 


Once a deviation of the Higgs couplings to the SM is observed, simultaneous measure- 
ments of various couplings can be used to distinguish different types of 2HDM. In Fig. 4, we 
plotted Ak, vs. A(ky/Ka) for four different types of 2HDM, for |cos(8 — a)| = 0.01 (green), 
0.02 (blue) and 0.03 (red). Left and right panels are for negative and positive cos(G — a), 
respectively. Dashed lines are for different values of tan 6, as labeled next to the lines. We 
also indicated the experimental precision of those couplings by the black cross, with a random 
central point, for the purpose of comparison. 

Type-I, II, L and Type-F are well separated, occupying the second, fourth, first and third 
quadrant, respectively, for cos(G — a) < 0. For the cos(G — a) > 0, the behavior is similar, 
except for exchange of the quadrants of Type-I + Type-II and Type-L + Type-F. For Type-I 
and II, Ake = Akg. Therefore, the variations of the coupling deviation with cos(G — a) are 
small, given the small allowed cos(8 — a) range that we picked. Smaller tan 8 value leads to 
larger deviation from the SM value for the Type-I, while the opposite occurs for Type-II. For 
the Type-L and Type-F, Ake and A(Ky/kq) spreads over the whole region, depending on the 
values of tan J and cos(8 — a). 


5.2 2HDM loop-level results in the alignment limit 


Other than the tree level deviation of the light Higgs couplings (as well as the loop generated 
hgg and hyy couplings with SM particles) in the 2HDM away from the alignment limit, 
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heavy Higgses in the 2HDM could also provide loop corrections to those couplings [68-70]. 
While these contributions are typically small given the loop-suppression and heavy Higgs 
masses suppression, they become the dominant correction to Higgs physics in or close to the 
alignment limit, cos(3 — a) ~ 0. Heavy Higgs contributions to hgg and hyy could also be 
important, given the loop suppressed SM values at the leading order. 

In this section, we analyzed the implication of Higgs precision measurements on the heavy 
Higgs loops, and explored the sensitivity to the heavy Higgs masses, as well as the Higgs self- 
couplings that enter the loop corrections. For simplicity, we set my = my+ = ma in the 


following discussion, which satisfies the EW precision p-parameter constraint automatically. 
We also worked under the most challenging scenario of 2HDM with tree-level alignment limit 
cos(( — a) = 0 to show the relevance and importance of these loop corrections. 


5.2.1 Theoretical and experimental constraints 


Heavy Higgs loop corrections would involve the Higgs self-couplings, which are constrained 
by various theories considerations and experimental measurements, such as vacuum stability, 
perturbativity and unitarity, as well as heavy flavor, electroweak precision measurements, 
and LHC direct searches. We briefly summarized below the constraints we adopted in our 
analyses. 


e Vacuum stability 


In order to have a stable vacuum, the following conditions on the quartic couplings need 
to be satisfied [71]: 


A, > 0, AQ > 0, 4X3 > — vy À1A2, Ag + à4 — As] > —vV/ à1A2. (5.8) 


e Perturbativity and unitarity 


We adopted a general perturbativity condition of |A;| < 47 and the tree-level unitarity 
of the scattering matrix in the 2HDM scalar sector [72]. 


EW precision measurements from the LEP constrain the p-parameter, or equivalently, 
the amount of custodial symmetry breaking in the 2HDM. Studies [73, 74] showed that the 
charged Higgs mass is constrained to be close to the mass of either of the neutral Higgses (H 
or A) in order to satisfy the EW precision measurements. In our analyses below, we adopted 
the simplification of my+ = my = ma = mg so that p-parameter constraint is automatically 
satisfied. The results we obtained below show characteristic features of the Higgs factory 
sensitivities to heavy Higgs mass, even though the numerical values might differ if masses for 
heavy Higgses deviate from the simplified relation. 

Under the assumption of alignment limit and equal mass for all the heavy Higgses, all 
the Higgs quartic couplings are related to a particular linear combination of més, and ii 


Av? = ms — m?,/(sin B cos 8) and we have Av? = \3v? — mẹ = —Aqv? = —A5v?. The above 
theoretical considerations can be translated to 
m? + Av? tan? 8 > 0 5 (5.9) 
E , h ' tan? B , 


=ņ17= 
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Figure 5. The shaded region indicates the surviving region of 2HDM parameter space of tan 8 
vs. Av?, after vacuum stability (blue lines), unitarity and perturbativity (red lines) are taken into 
account. The corresponding values of A are shown in the upper axis. Here we took the simplification 
of my+ = My = ma = mọ and alignment limit of cos(3 — a) = 0. 


dv?(GeV2) | —1002 | 0 | 100? | 300? | 500? 
tan Bmin 0.80 0 | 0.14 | 0.43 | 0.73 
tan Bmax 1.25 | +00 | 7.01 | 2.33 | 1.37 


Table 4. A few benchmark values for \v?, and the corresponding range of tan 6 given the theoretical 
considerations. 


for Av? < 0, and 


1 64r?vt + 5m} — 48ru? m? + 8A204 — 4m? dv? 


tan? < , 
aes tan? 8 3Av? (872 — 3m?) 


(5.10) 


for Av? > 0. In Fig. 5, we show the allowed shaded region in tan 8 vs. Av? plane given the 
theoretical considerations. Region in tan 8 and 1/ tan 8 is symmetric, which is obvious from 
Eqs. (5.9) and (5.10). A few representative values of Av? that we used in our later analyses 
and the corresponding acceptable region of tan 8 are shown in Tab. 4 as well. Note that 
for A = 0, i.e. m3 = m?,/(sin 8 cos 8), tan 8 is unconstrained, which is consistent with the 
results of Ref. [73]. Given the symmetry between tan 8 and 1/tan p in the above conditions, 
the largest region on Av? occurs at tan 8 = 1: 


— m? < Mv’ < (600 GeV)’, (5.11) 


which gives —0.258 < ÀA = —A4 = —A5 < 5.949 and 0 < A3 < 6.207. 
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There are direct searches of non-SM heavy Higgses at the LHC [75], with the dominant 
search channel being A/H — rr. In the framework of MSSM, m,/# is excluded to about 
250 GeV for tan 8 > 1.0, and about 1.5 TeV for tan 8 > 45. Since the branching ratio of 
the dominant search channel A/H — rr could be highly suppressed once other exotic decay 
channels of the non-SM Higgs opens up [76-78], the current exclusion limits could depend 
highly on the non-SM Higgs spectrum. The direct search limits on the heavy charged Higgs 


H= above m; are relatively weak given the large SM backgrounds for the dominant H= — tb 


channel, and the relatively small branching fraction of H* — rv [79]. 

Flavor physics consideration usually constrains the charged Higgs mass to be larger than 
about 600 GeV for the Type-I] 2HDM [79]. However, the charge Higgs contributions to 
various flavor observables can be cancelled by other new particles in a specific model [80] 
and be relaxed consequently. In our analyses, we focused on the indirect search potential 
of the Higgs factories on the masses of heavy Higgses. Therefore, we did not impose flavor 
constrains on the 2HDM parameter space, as well as the LHC direct search limits. 


5.2.2 2HDM loop effects 
We defined the normalized Higgs coupling including loop effects as: 


2HDM 2HDM 
2HDM _ Itree + Joop 


loop SM SM °? 
Gtree pg Noop 


(5.12) 


2HDM 


in which gop is 2HDM loop correction involving both SM loop corrections and non-SM 
parts. 


To the leading order of 1-loop correction, Eq. (5.12) simplifies to 


2HDM _ „tree 2HDM 
K1—-loop — * F AK1 loop: (5.13) 
: tree — „2HDM /„SM 2HDM — „2HDM /„SM 2HDM ; : 
with K = gee /Gtreer AKI loop = J1—loop/ Itree Where g] loop is the one loop corrections 


with non-SM particles running in the loop. In the alignment limit of cos(3—a) = 0, Kte = 1, 
the expression is simplified to be 


2HDM _ 2HDM 
K1- loop lalignment =1+ Aki bop’ (5.14) 


The expressions of the non-SM Higgs loop correction to the Higgs couplings are summa- 
rized in Appendix B [68-70] under the alignment limit cos(8 — a) = 0 and with the mass 


simplification relation my+ = Mmao = myo = mg. Note that the tree level relations of 
Kw = kz and Kk, = K, are still approximately valid at 1-loop level. 

Following the same global fitting techniques as in the tree level case, we obtained the 
95% C.L. constraints on the tan 8 vs. mg plane for four different types of 2HDM given CEPC 
precisions, as shown in Fig. 6. Four benchmark values of VAv2 are chosen: 0, 100, 300, and 
500 GeV, which correspond to red, blue, green, and orange curves. Regions to the right of 
the curves for large mg or to the left of the curves for small mg (enclosed region for the blue 
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Figure 6. 95% C.L. constraints on tan 6 vs. mg plane based on CEPC Higgs precision measurements. 
The orange, green, blue and red (from right to the left at the large mg region) are for VAv? = 


2 mip 
MG ~ sacp 
to the left of the curves for small mg (enclosed region for the blue curves) are allowed. Four panels 


are for Type-I, II, L, and F 2HDM, as labeled at the top of each figure. 


=500, 300, 100, 0 GeV, respectively. Regions to the right of the curves for large mg or 


curves) are allowed by the CEPC precision measurements assuming no deviations from the 
SM predictions are observed. 

For Type-I, small values of tan 8 S 0.5 are excluded since all the non-SM Higgs Yukawa 
couplings are proportional to 1/tan(. The dependence on tan 8 is weak once tan = 2. 
While for smaller values of \v?, Mg as low as 125 GeV are allowed, mg is constrained to be 
larger than 500 GeV for VAv2 = 300 GeV, and 1200 GeV for VAv2 = 500 GeV, given that 
heavy Higgs loop corrections enhance as Av? increases. The small kink around Mg ~ 350 
GeV is due to the top quark threshold effects in Yukawa couplings. The relaxed constraints 
on tan f around mg ~ 500 GeV are caused by the smallness of Akb cr, which flips sign near 
that region. 

For Type-II, while the generic features of the excluded region are the same as those in 
Type-I, there are three major differences. Firstly, the constraints at large tan 8 get tighter 
since both the down-type Yukawa and lepton Yukawa (in particular, those of bottom and 
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tau) are proportional to tan 8. Therefore, loop contributions are more tightly constrained at 
large tan 8. Secondly, the top quark threshold effects are stronger since the relevant terms 
do not have the 1/tan? 3 suppression as in Type-I. Thirdly, the constraints at small tan 8 get 
weaker. This is because only up-type Yukawa couplings are proportional to 1/ tan 8, while 
the precision on that is worse than the bb and r7 channels. 

The results for Type-L are similar to those of Type-IJ, with small tan 8 constraints 
getting stronger since the down-type Yukawa couplings are now proportional to 1/ tan 8 as 
well. The sensitivity at large tan 8 becomes a bit weaker since the lepton Yukawa coupling 
structure is the same as that of Type-II. The results for Type-F are almost the same as that 
of Type-II, since for the most precisely measured couplings, hZZ and hbb, the dominant loop 
contributions (from bottom quark and heavy Higgses) in these two types are identical. 
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Figure 7. 95% C.L. Higgs precision constraints on tan 6 vs. mg plane, assuming integrated luminosity 
of 2.5, 5, 10, and 25 ab~’. Results for Type-I (top panels) and Type-II (bottom panels) are presented, 
for two different benchmarks of v Av? = 0 (left panels) and 300 GeV (right panels). 


To explore the effects of different running scenarios with different integrated luminosity, 
we presented the 95% C.L. Higgs precision constraints on tan 8 vs. mg plane in Fig. 7, with 
a rescaled error bar of different integrated luminosity of 2.5, 5, 10, and 25 ab~!, based on 
CEPC 5 ab”! results. Here we have assumed that statistical error dominates[1, 39]. Both 


A 


the upper and the lower limits on tan 8, as well the lower limits for mg at large mg region 
get stronger with increasing luminosity. For both Type-I and II with VAv? = 300 GeV, the 
lower limit on mg increases from 400 GeV to 700 GeV when the luminosity increases from 
2.5 to 25 ab 
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Figure 8. 95% C.L. Higgs precision constrained regions in mg vs. tan for CEPC (red lines), ILC 
(blue lines), and FCC-ee (green lines). Left panel is for Type-I and right panel is for Type-II with 
Vv? = 300 GeV. The CEPC results using «-fit are shown in red dashed lines. 


To compare the sensitivity of different future Higgs factories with the running scenario 
listed in Table. 1, as well as the effects of Au/u-fit vs. «-fit, we show the 95% C.L. constrained 
region in tan 8 vs. mg plane for CEPC (red lines), ILC (blue lines), and FCC-ee (green lines). 
For CEPC precision, the Ay/,u-fit results are shown in solid lines and «-fit results are shown in 
dashed lines. ILC and FCC-ee have similar sensitivities, both better than the CEPC results. 
In addition, results with Ay/ fit are better than the results with «-fit, confirming the tree 
level fitting results that including the correlations between different couplings is important. 
Overlooking those correlation effects might lead to overly conservative results. 

Results we presented so far are based on the assumption of a fixed value of Av? = 
més, — m7,/(sin8 cos 8). It is convenient to do so since Av? directly enters the tri-Higgs 
couplings. Another common strategy is to fix the value of the soft Z2 breaking parameter 
m?,. The global fitting results are shown in Fig. 9 for fixed values of m7, to be 0 (red), 150? 
(blue), 500? (green) and 1000? (orange) GeV”. As before, the top quark threshold effects 
show up in the blue region. Negative values of m?, typically lead to large Av’, therefore 
much worse values of x? and disfavored. The survival regions now exhibit (almost) enclosed 
behavior, with both upper and lower limit on mg for a given value of tan 3. The location of the 
allowed mg region also shifts to the larger mass for larger values of m12. For the asymptotic 
large mg and tan @ region, the allowed region centered around ms = m?,/(sin 8 cos 8) line, 
which minimizes the corresponding v?. 

To study the sensitivity of Higgs precision measurements on the model parameters, as 
well as to demonstrate how well Higgs precision measurements can be used to distinguish 
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Figure 9. The same global fitting results with CEPC as Fig. 6, with different values of mî., instead 
of \v?. The survival regions exhibit (almost) enclosed behavior, with both upper and lower limit on 
mg for a given value of tan £. 


different types of models once certain deviations in the couplings are measured, we show the 
loop corrections of Ax, vs Aw for four different types of 2HDM, for mg = 600 GeV, with 
Vv \v2 = 300 GeV (left panels), and 500 GeV (right panels). Different values of tan 8 are also 
indicated in the curves with dots to show the sensitivity of Kb,- to the values of tan £. 

For Type-I, the loop corrections to both «Kp and «, are large for small tan 6, about 4.3% in 
Ak, and 2.3% in Ax, for tan 8 =0.7 with Vv? = 500 GeV, to be well explored by precision 
measurements. The contribution quickly decreases once tan 8 gets larger. For Type-II, the 
dependence on tan is flipped, with the corrections to «, significantly larger than that of 
Type-I. For Type-L, the correction to «p is large at small tan 8, while the corrections to kK, 
increases for larger tan 8. The opposite behavior occurs for Type-F. For larger values of V Xv , 
the loop corrections get even larger. 

Since each type of 2HDM swipes out a different region in Ax, vs. Ax, space, simultaneous 
measurements on Kp and &, with the estimated precision could help to distinguish four types 
of 2HDM models. For CEPC with 5 ab~! 7-parameter fit, the estimated precision is about 
1.2% for k, and about 1.3% for «,, as indicated by the cross in Fig. 10, which provide enough 
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Figure 10. Loop corrections of Ax, vs Ax, for four different types of 2HDM, for my = 600 GeV, 
with v Av? = 300 GeV (left panels), and 500 GeV (right panels). The values of tan 8 are indicated in 
the plots with dots. The black cross indicate the estimated experimental errors with a random central 
point. 


precision to separate contributions from different models. Note that without knowing the 
value of tan 8 a priori, there is degeneracy between Type-L and Type-F, in particular, if 
only deviations in ky, is observed while «x, deviation is small. However, with the theoretical 
constraints discussed in Sec. 5.2 imposed, a constrained range of tan 8 as shown in Fig. 5 
and Table 4 could help to lift the degeneracy. Since the loop corrections typically get smaller 
for larger values of mg and smaller values of A, the potential of kp and k, measurements to 
determine the values of tan 6 becomes limited accordingly. 


6 Composite Higgs models 


The class of BSM physics that features strong dynamics genuinely leads to testable predica- 
tions through precision measurements. Composite Higgs Models are very representative class 
of models that ties directly to the hierarchy problem. The Higgs boson, instead of being a 
fundamental particle as in the SM, SUSY model or gauge extensions of the SM, would be a 
composite particle as a bound state of additional strong dynamics of the underlying physics 
model. Higgs now can be viewed as a pseudo-Nambu-Goldstone boson, whose mass is pro- 
tected by the global symmetry breaking scale parameter f generated by the condensation of 
strongly interacting particles. The separation between the electroweak scale and the com- 
posite scale f is naturally a tuning parameter, representing the fine-tuning of the underlying 
composite model. 

There is a vast range of plausible composite Higgs models. To evaluate the physics reach 
of the Higgs precision measurements on this broad class of models, we adopted two approaches 
in discussing the results. The first approach is to interpret the Higgs precision results in the 
Minimal Composite Higgs Model (MCHM) [81] with various embedding of the partners of 
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the SM fermions’. The second approach is to adopt the language of EFT and followed the 
discussion of the patterns of strong interactions [84], comparing the implication of the Higgs 
precisions on the so called Accidentally Light Higgs (ALH), Strongly Interacting Light Higgs 
(SILH) [85], and a general SILH (GSILH). Each of those patterns has different assumptions 
on the symmetries of the underlying strong dynamics. We discussed in details about their 
underlying assumptions in the following sections. 


6.1 Minimal Composite Higgs Models 


The MCHMs represent a minimal embedding of the Higgs as a pseudo-Nambu-Goldstone 
boson of the global symmetry SO(5)/SO(4) that respects custodial symmetry. We investi- 
gated and found a way to best present limits on various tower fermion representations in the 
MCHMs, following the notation and calculations detailed in Ref. [86]. 

In the minimal coset SO(5)/SO(4), gauge invariance fixes the rescaled gauge coupling 
Ky to be 


g% 
ky = a =y1-£, (6.1) 


hVV 
where £ = v?/ f? parameterizes the vacuum misalignment. To simplify the notation, in this 
section, we referred to the rescaled couplings of the SM Higgs Ki as ki. The modification to 
the Yukawa couplings depends on the fermion representations. In many models, the rescaled 
Higgs to fermion couplings k+ or Kp, can be either 


ES =, Fe=4/1—€, (6.2) 


if summed over all the tower contributions. In these cases € is the only model parameter at 


leading order for a specific model. For some more complex embedding of the fermions, such 
as MCHMy4- 14-190, MCHM 4-519 and MCHM3_ 44-10, the Higgs couplings, are modified in 
a more complex way [86]. The corresponding Yukawa couplings are related to functions F345, 
which depend on several microscopic Yukawa couplings, with the explicit expressions can be 
found in Ref. [86]. The resulting Higgs couplings vary in a certain range, even for a fixed 
value of €. In several limiting cases when one of the microscopic Yukawa couplings turns off, 
these additional coupling functions reduce to simpler functions of Fy and F>. 

Note that these simple closed form expressions are obtained by summing over the infinite 
tower fermions. In reality, the summation is truncated after a few tower fermions, generating 
a more complicated and scattered relation between model parameters. For simplicity and as 
the modifications to the Higgs couplings are dominated by the first few tower fermions in 
most cases, we adopted these simplified formulae to obtain a qualitative physics reach of the 
Higgs precision program in the MCHMs. 

We tabulated in Table 5 various MCHM fermion representations and the corresponding 
leading modifications to Ky (so as to Kg) and Kp, while the modifications to Kz always follow 


>For recent studies focus on the future collider perspective, see, e.g., Refs. [82, 83]. 
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5, 10 
MCHM | 14-1-10 5-10-10 
Reva: 14-10-10 10-5-10 | 5-5-10 5-1-10 14-14-10 14-5-10 5-14-10 
10-14-10 
EA F, F, Fi B T; F, F; 
Kp F, Fi Fə Fy F, Fi PF, 
CEPC 
Ex 10° 2.56 2.36 4.19 3.87 2.78 — 2.56 | 2.71 — 2.36 | 2.36 — 2.04 
f [TeV] 4.86 5.06 3.80 3.95 4.67 — 4.86 | 4.72 — 5.07 | 5.07 — 5.45 
ILC 
Ex 10° 2.19 2.02 3.44 3.20 2.31 — 2.19 | 2.06 — 2.01 | 1.87 — 1.72 
f [TeV] 5.26 5.48 4.19 4.35 5.12 — 5.26 | 5.42 — 5.48 | 5.69 — 5.93 
FCC-ee 
¿x 10° 1.80 1.66 3.06 2.74 1.85 — 1.80 | 1.70 — 1.66 | 1.66 — 1.41 
f [TeV] 5.79 6.04 4.45 4.70 5.72 — 5.80 | 5.97 — 6.05 | 6.05 — 6.56 


Table 5. 95% C.L. bound on € and global symmetry breaking scale parameter f for the MCHMs 
with various embedding of the fermion content, under CEPC, ILC, and FCC-ee Higgs precisions. 


function Ff. All the Higgs couplings only depend on one parameter € for the first four cases 
when «k and «p is either F\ or Fy. The remaining three embeddings have extra parameter 
dependence that enters F345. 

Adopting the global fitting method as described in Sec. 3, we obtained the 95% C.L. 
range on the €, which is listed in Table 5 for CEPC, ILC, and FCC-ee Higgs measurement 
precisions. We further translated the corresponding limits on € into the more intuitive limits 
on composite scale parameter f. We note that «į can be well constrained by h —> gg here 
as we assume the h — gg process does not receive additional new physics contribution. We 
did not make any assumption on the model prediction of ke. If further assumptions are made 
(e.g. Ke = Kt), a Marginal improvement on the constraint of € is expected, as the bounds 
on Ky and «p are much more constraining. At 95% C.L., € can be limited to be less than 
a few times 1073, assuming future measurements follow SM expectations. The composite 
scale is constrained to be bigger than 4 TeV. For the embeddings with «+ related to function 
F345, the 95% C.L. bounds of € and f vary in a certain range, given the extra parameter 
dependence. 

Once deviations of the Higgs couplings are observed at future colliders, different embed- 
ding of the fermion tower might be distinguished through the predicted correlations between 
Higgs couplings. In Fig. 11, we show in the Aky-Ak; plane the 68% C.L. and 95% C.L. 
exclusion limits of the Higgs precision measurements of CEPC, FCC-ee and ILC, obtained 
from a global fit in a reduced parameter space of the four parameters Ky, Kt, Ke, and Ky. The 
loop contributions of these parameters in the hgg, hyy, hZy couplings are also included. In 
the right panel, the 95% C.L. from the LHC 300fb~! and 3000 fb~! runs are also shown for 
comparisons with the lepton colliders. 
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Figure 11. The 68% C.L. (solid lines) and 95% C.L. (dashed lines) constraints in the (Aky, Ars) 
plane at various future Higgs factories from a four-parameter fit with ky (= Kz = Kw), Ke, Ke, and 
Kp. The left panel shows the constraints for CEPC, FCC-ee and ILC, while the right panel shows the 
LHC 95% C.L. constraints as well. Also shown are the predicted deviation of k; and «y for different 
fermion embeddings, as a function of €. The magenta and cyan lines correspond to kK; = F and Fo, 
respectively, while x = F3,4,5 are covered by the gray region. The labels on each part of the gray 


region indicate which ones of Kk, = F345 cover this part. 


The four-parameter fit is very useful in capturing the main characterization of the MCHMs 
without making specific assumptions on the fermion representations. The projection to the 
Aky-Ak; plane allows us to focus on the two couplings most relevant for the MCHMs. For 
specific fermion representations, the parameter space is further constrained, implying certain 
correlations between Axy and Akg. For kų = Fi or Fo defined in Eq. (6.2), both Akz and Akg 
are fixed by the value of €, shown by the magenta and cyan lines in Fig. 11. For Ky = F345 
with the variation of additional model parameters, the predicted range is covered by the gray 


region. 


The right panel of Fig. 11 shows clearly that the Higgs measurements at future lepton 
colliders can significantly improve the constraints on the MCHMs from the corresponding 
measurements at the LHC for more than one order of magnitude. Once deviations of Higgs 
couplings are measured, different embeddings of Fermion contents could also be tested. This 
is, of course, due to the much better determination of the Higgs couplings at lepton col- 
liders. In particular, the coupling of Higgs to the Z boson is very well constrained by the 
Higgsstrahlung processes, while at the LHC it is probed by the Higgs decay h — 4£ or the 
VH production mode and suffers from systematic uncertainties. 


— 927 — 


Ox = 5(0,|H? |)" Oca = gH GA GAM 
Ow = 8(Htot D!H)DYWS, | Oy, = YalH|? Qr Hur 
Op = Z (HYD" H) By, Oy, = YA Ordr 

Oyuw = ig(D"H)to°(D”H)W}, | Oy, = Ye|H|?LrHer 

Onp = ig'(D“H)'(D’ AH) By Osw = 79€aveW a” We, Were 
Ose =o" |b? Bye’ 


Table 6. Operators in the SILH basis considered in our study. Assuming no corrections to electroweak 
precision observables, only one combination of Ow and Og is left after imposing the electroweak 
precision constraints of cw +cg =0. While Y,,/a/e are 3 x 3 matrices in general, we only consider the 
relevant diagonal ones Yt, Yc, Yb, Yr and y,,. In this table, egy. is the totally anti-symmetric tensor for 
SU(2). 


6.2 General Patterns of Strongly Interacting Light Higgs 


Given the richness and depth of strong dynamics with a light Higgs, it would be very infor- 
mative to understand Higgs factory potentials on strong dynamics models under some generic 
arguments about the patterns of corresponding EFTs. In our analyses, we considered three 
typical patterns of strong dynamics with a light Higgs [84]: ALH, SILH and GSILH. These 
three cases are categorizations of general patterns of strong interactions rather than explicit 
models. 

To capture the main feature of these cases without the loss of generality, we worked in 
the EFT framework with dimension-6 (D6) operators, parameterized by 


1 
Le = m2 D ciO; ; (6.3) 


where m, is the new physics scale and the coefficients c;s are generally determined by the 
coupling of new physics, denoted as gx, as well as the SM gauge couplings (gs, g and g’, for 
SU(3), SU(2) and U(1), respectively) and Yukawa couplings yy. A collection of the relevant 
operators ©; are listed in Table 6. Amongst these operators, Oy requires renormalization 
of the Higgs field and shifts the Higgs couplings universally to other SM particles, Oy oper- 
ators further modifies the Higgs Yukawa couplings. Operators Ogg, OHB, Onw, and OGG 
directly modify Higgs to gauge boson couplings. Operators Ow, Op, Onw, OHB, and Osw 
also contribute to electro-weak precision observables and anomalous TGCs. Consequently, 
having additional measurements other than Higgs properties is crucial for developing global 
constraints on these operators. 

For three strong dynamics cases that we considered here, the estimated parametric count- 
ing and scaling of the operators are listed in Table 2 of Ref. [84], which are reproduced here 
in Table 7. For each of these cases, there are only two free parameters: g, and m,. However, 
it should be noted that Table 7 provides an estimation on the size of the c;’s rather than 
their predicted values, which are available only if the model is specified. We therefore do not 
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On | Ow | Op | Onw | Ous | Ops | Oae | Oy, | Oya | Oye | Osw 
ALH | 2/1 [1 1 1 1 rAr 
Z 2 5 
GSILH | g | 1 | 1 1 t geet |aeee| | Ge |e = 
7 z 7 7 $ 

SILH | g | 1 1 | e | se | we | wr | 9 | 9% | 9% | te 


Table 7. Estimation of the magnitudes of c;’s for the operators in Table 6 for the three scenarios 
of ALH, GSILH, and SILH. Note the difference between this and Table 2 of Ref. [84] is due to the 
different normalization of operators. 


assume any relations among the operators in Table 7, but only use the estimated coupling 
size to derive the scale of new physics mą. 

We derived the reach of the scale of new physics m, as a function of g, for three individual 
scenarios from the 95% C.L. constraints on the overall coefficient of EFT operators derived in 
Ref. [34], which was translated to m,/,/cj as shown in Fig. 12 for CEPC (red) , ILC (blue), 
and FCC-ee (green) Higgs precisions. The Higgs measurements from HL-LHC are included in 
the total y? to further optimize the reach. The HL-LHC improves the limits through better 
diphoton statistics and ttH processes that compliment Higgs factory measurements. Note 
that this is different from the global analyses on the 2HDM and MCHM in precious sections, 
where the inclusion of HL-LHC Higgs measurements does not affect the results much. In 
addition to the Higgs measurement inputs as listed in Section. 2, the angular observables in 
ete — hZ and the constraints on anomalous TGCs from measurements of e'e~ — WW 
are also included in Ref. [34] to help discriminate different EFT parameters and maximize the 
overall precision reach. Ref. [34] also assumed that CEPC is able to collect 200fb~! data at 
350 GeV. Both CEPC and FCC-ee 350 GeV measurements are obtained by scaling from ILC 
350 GeV, assuming statistical uncertainties dominate. The difference between unpolarized 
beams and polarized beams has also been taken into count. 


Note that Ref. [34] only focuses on operators in Table 6 that can be probed by Higgs 
processes and diboson production (ete — WW). Furthermore, relation cw + cg = 0 is 
imposed such that there is no additional contributions to the electroweak precision observ- 
ables, resulting in only one combination of Ow and Og survives. Ignoring the flavor changing 
effects and only taken into account the relevant diagonal Yukawa couplings Yt, Ye, Yo, yr and 
Yu, we ended up with 12 coefficients to be constrained by the global fit. 


Fig. 12 shows the results from individual constraints on operators obtained by switching 
on one of them at a time (light shade) and the other from a 12-parameter global fit (solid 
shade). The single operator fit results are typically factors of a few better for operators 
involving electroweak gauge bosons, due to their correlation in Higgs physics. We also note 
here the typical scale from a global fit on these operators are around 1—3 TeV with cœ ~ 1, 
with the exception of m,/cgp and m,/cgqg, which could reach around 10 TeV. This is because 
the leading order contribution from the SM to h > yy, gg appears at one-loop level. 


Using the bounds on m,/,/c; in Fig. 12, the reach of the new physics scale m, can be 
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Figure 12. The 95% C.L. constraints on the overall coefficient of O; from Ref. [34], which are 
translated into the SILH basis and presented in the form m,/,/c;. Estimated Higgs measurement 
precision from CEPC, ILC, and FCC-ee are used with the inclusion of HL-LHC Higgs precision. 
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Figure 13. The 95% C.L. limit on the new physics scale for the three cases of ALH (dotted lines), 
GSILH (dashed lines) and SILH (solid lines) as a function of g, with CEPC (red), ILC (blue) and 
FCC-ee (green) precision. The left plot is the results obtained using fit on individual operators, and 
the right plot is the results obtained using global fit of all 12 operators. The operator that is most 
sensitive (therefore determines the best reach) is labelled alongside the curves. 


derived as a function of the strong coupling gx. The reaches are presented in Fig. 13 at 95% 
C.L. for three cases under future Higgs factories. We varied g, from 1 to 10, to cover the 
typical range of strong interaction couplings 1 < gą < 4r. Two sets of bounds are shown 
for each of the three cases of ALH, GSILH and SILH, one from individual fit (left panel) 
and the other from a global fit (right panel). The constraints obtained from individual fit 
provide a optimistic estimation as in realistic models it is unlikely that only one D6 operator 
is generated. The global fit, on the other hand, provides a conservative bound, with all 
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coefficients of the relevant D6 operators treated as independent free parameters. For each 
cases, the operator that provides the strongest constraints is labelled alongside the curves. 

For gą varies between 1 to 10, the reach in m, is about 5 to 50 TeV for individual fit, 
and about 2 to 30 TeV for global fit. For individual fit, the CEPC and ILC reach is similar, 
except for small g, region, while FCC-ee reach is slightly better. For global fit, the ILC reach 
improves over CEPC and FCC-ee, due to the additional measurements of Higgs and diboson 
processes at higher center of mass energy. While the reaches in GSILH and SILH are almost 
identical, the reach in ALH is quite different, in particular, on the Ogg operator where the 
loop suppression from the shift-symmetry argument is absent. 

The dominating Higgs production process, the Higgsstrahlung process ete” — hZ, is 
ideal for probing Oz, which provides a universal shift to all the Higgs couplings. However, in 
a global fit its bound suffers from a large degeneracy with the other parameters, in particular 
with operators that could generate anomalous hZZ couplings with Lorentz structures different 
from the SM one. Therefore, in individual fits Oy provides the best reach for GSILH and 
SILH, while for global fit, the best reach comes from O,, instead at CEPC and FCC-ee. ILC 
with the run at 500 GeV is more powerful at discriminating these different operators through 
the WW-fusion production of the Higgs boson and ZH associated production [87], therefore 
with a global fit Oy still has the best reach. For smaller values of gą, Ow provides the best 
reach since its coefficient is independent of gx. 

For ALH, Ogg provides the best sensitivity to m,, around 50 TeV in the individual 
operator fit case, and around 20 TeV in the global fit case for gx <6. This is because OGG 
can be generated at tree level in ALH, while at one loop level in the SM, GSILH and SILH. 
In the global fit, the bound on Oga usually suffers from a flat directions associated with 
Oy, which also contributes to the hgg coupling by modifying the top Yukawa coupling in the 
loop. Therefore, in ALH, the most sensitive operator in global fitting at large gx is Oy, at 
CEPC/FCC-ee and Oy operator at ILC, similar to the case of GSILH and SILH. 

With these results on the typical patterns of strongly interacting light Higgs, we can see 
how different components of the Higgs precision programs compliment each other in probing 
models of strong dynamics. In particular, running at different center of mass energies of Higgs 
factories can be sensitive to a particular set of operators. 


7 Conclusion and outlook 


In this work, we studied the implication of the Higgs precision measurements from various 
Higgs factories (CEPC, ILC, FCC-ee) on new physics models. In particular, we focused 
on two types of new physics scenarios: perturbative models, and strong dynamics. For 
perturbative models, we studied the SM plus a singlet extension as well as the 2HDM as 
explicit examples, and explored both the tree level effects as well as loop corrections. For 
strong dynamics, we studied the MCHM as a concrete example, and also explored the EFT 
language in three typical patterns of operators that arose from strong interactions. We took 
the estimated precision from Higgs factories, as well as LHC results, and performed a global 
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fit in the relevant parameter space to derive the reach in different models. In the case when no 
deviation from the SM predictions is observed, a 95% C.L. reach in the parameter range was 
derived. In the case when deviations from SM predictions are observed, we further explored 
how different models can be distinguished. 

For the SM with a real scalar singlet extension, we obtained the 95% C.L. limit of the 
singlet-doublet mixing angle sin 0 to be 0.062, 0.058 and 0.052 for the CEPC, ILC and FCC- 
ee, respectively, at 95% C.L. We also studied the more general case with the induced Oy and 
Og operators, and obtained constraints on the Wilson coefficients cy/A? and c¢/A?. While 
cy /A? is tightly constrained to be around 10~?/v?, the limit on cg/A? is about an order of 
magnitude weaker. 

For the 2HDM, we analyzed four different types of 2HDM: Type-I, II, L and F. For 
tree level effects, we found that |cos(8 — a)| is tightly constrained to be about 1074 or 
smaller at small and/or large tan 8 for Type-I, L, F, and reached the maximum range of 
about 0.005 — 0.007 for tan ~ 1. For Type-I, large tan 6 region is much less constrained: 
|cos(3 — a)| < 0.08. For 2HDM loop effects under the alignment limit, the lower bounds on 
the heavy Higgs masses depend on the Av? = m3 — mî, /(sgcg): about 350, 450, and 1100 
GeV for V\v? = 100, 300, 500 GeV, respectively. For smaller values of V v2 < 100 GeV, 
a small mass region of 125 GeV < mg < 350 GeV is viable as well. We also explored the 
constrained parameter space in tan vs. mg plane when mî is varied as an independent 
parameter. 

For the MCHM, the main parameter is € = v?/f?, which is linked to global symmetry 
breaking scale f. With different embeddings of the fermion content, we obtained the 95% 
C.L. limits on € to be about a few x1073, corresponding to f > 4 TeV. We also studied the 
ALH, GSILH and SILH as three representative patterns of EFT operators in strong dynamics. 
The 95% C.L. constraints strong dynamics scale m, as a function of strong coupling g, were 
obtained, either with an individual fit to a single operator, or through a global fit of all twelve 
operators. For gą varying between 1 to 10, m, varies between 5 to 50 TeV for individual fit, 
and 2 to 30 TeV for global fit. 

We also studied the reach for three different Higgs factory machines: the CEPC, ILC and 
FCC-ee under typical running scenarios as summarized in Table 1. The reaches of FCC-ee 
and ILC are slightly better than that of CEPC. With higher center of mass energy running 
of ILC, it has advantage of being sensitive to operators that contribute to hWW couplings. 

In our analyses, we used the global fit to the signal strength u; for each individual search 
channel. This has been shown to provide stronger constraints on parameter spaces over the 
usual «-scheme when correlations between «is are not treated properly. 

To summarize, the impressive precisions in the Higgs measurements that can be achieved 
in future Higgs factories provide us a powerful tool to explore the new physics models through 
indirect probe. Those indirect searches to new physics are complementary to direct searches 
at current and future colliders, and sometime even supersede direct limits given the hadronic 
environment at the LHC and the limited energy reach of lepton colliders. Our studies pro- 
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vided concrete examples of how sensitive Higgs precision measurements can be to particular 
model parameters, new particle masses, or symmetry breaking scales. While our results are 
specific towards the SM singlet extension, 2HDM, MCHM, and three particular types of 
strong dynamics set up, the approaches we took are more general and can be applied to any 
model set up, as long as Higgs couplings and other relevant interactions are specified. Our 
study also provides guidelines to future Higgs factory design, including choices of center of 
mass energies and the corresponding luminosities. 
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A LHC Higgs measurements 


In Table 8, we listed the LHC Run-I Higgs measurements for yy, WW, TT, bb and ZZ for 
various production channels. In Table 9, we listed the estimated precision on normalized 
Higgs signal strength for 300 and 300 fb~‘ integrated luminosity, as well as the corresponding 
production contamination. These two tables are used as input for our global fit using LHC 
data. 


B 2HDM loop corrections 


The formulas of the one-loop correction from the 2HDM to the Higgs couplings are summa- 
rized here, under the alignment limit cos(G — a) = 0 and with the universal mass assumption 
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Channel | Production | Run-I | Channel | Production | Run-I 
yy ggh 1.104933 TETT ggh 1.0706 
¥0.5 +0.4 

VBF 1.37555 VBF 1.3104 

+13 +14 

Wh 0.57 )'5 Wh —1.4°7%4 

F3. ¥2.2 

Zh 0.5439 Zh 2I 

7 +16 7 ¥3.7 

tth 2.2113 E tth —1.913:3 

ww* ggh 0.841917 bb Wh 1:0193 
VBF tan Zh Oates 
Wh Lela. tth 15 as 
Zh 5.9728 Le ggh Mat 
P99 490.31 

tth 5.0113 VBF Ora, 


Table 8. The Higgs signal strength in various channels measured at the LHC Run-I [47]. Correlations 
between different Run-I measurements quoted in Fig. 27 of Ref. [47] are taken into account. 


MH = mM Ao = MHo = Mg. 
2uA, H+H- 1 
2HDM 2 hH+H 2 2.9 2 2 
Akiy,1-loop = —— 9 [1 + 2m4+C0(0, 0, mp, Myt: Ms, M+) F— 
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2 l Ly miz/ (seca) 1 (23 gia/(s8¢8) ) mi, 
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1 mīa/(sgcg) mi J] 1 
73—9 ; B.1 
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with Abn = 6.53, representing the one-loop induced hyy coupling in the SM, which includes 
the W boson, top quark and bottom quark contributions. 


2M2, 
a cot 8 [m — 2M4)C(M4, m4, me; Mp, Meo, Mfr) 


1 
2HDM n 
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: Au/u True Origin 
Decay | Production 

300 fb7! | 3000fb-! | ggF VBF | WH | ZH | ttH 
comb 0.09 0.04 aan ee IE Eo 

Oj 0.12 0.05 am a 0. 0. 0. 

1j 0.14 0.23 im i 0. 0. 0. 

yy VBF-like 0.43 0.15 0.3(0.4) | 0.7(0.6) | 0 0 0. 
WH-like 0.48 0.17 0. 0. z 0. 2 

ZH-like 0.85 0.28 0. 0. i | Ss | se 

ttH-like 0.36 0.12 0. 0. so | om | # 
comb 0.07 0.04 Soe ca AE Seo 

VH-like 0.34 0.12 25 765 75 75 7055 
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VBF-like 0.33 0.16 ori ari at | mi Ta 
ggF-like 0.07 0.04 3309 3809 309 | 3809 | 3800 
comb 0.08 0.05 a 608 AE | ces 

N 0j 0.09 0.05 TE TA 0. 0. 0. 
1j 0.18 0.10 w o 0. 0. 0. 

VBF-like | 0.20 0.09 a 5o 0. 0. 0. 

Zy incl. 0.44 0.27 n pai 0. 0. 0. 
comb 0.26 0.12 E sia sooo. |) sean Seas 

bb WH-like 0.56 0.36 1 0. 0. 0. 0 
ZH-like 0.29 0.13 0. 0. 0 | @ | & 

tT | VBE-like 0.18 0.15 soa ae 0 0. 0. 
comb 0.38 0.12 T eo aa a ER 
PE incl. 0.45 0.14 172 1735 mE | is | T3 

ttH-like 0.72 0.23 0. 0. 0. 0. 1 


Table 9. The normalized Higgs signal strength uncertainty [48] and their corresponding production 
contamination information in various channels [88-94] measured at the LHC 300 fb~! and 3000 fb™t. 
Here the contamination information are for 3000 fb~', and almost the same for 300 fb~! except for 
the VBF-like H — yy process, which are shown in brackets. 


where f’ is the fermion whose electromagnetic charge is different by one unit from f, and 
Iş = +1/2 (—1/2) for f = u (d,e). 
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The scalar three-point couplings are given by 


2 2 
(m? + 2m3 = 272) TAD a m3 
5 ee. Agoo = 2 ot 28 (B.4) 
v v 


Àh = 


All the B and C function are Passrina-Veltaman Functions [95]. 
We should emphasize that all the above equations are in the alignment limit cos(S—a) = 0 


and mg = My+ = M40 = Myo. Following the fitting techniques at tree-level, kw = Kz and 


Ku = Kr Still hold at 1-loop. 
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